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Abstract 

We suggest a better mathematical method, fractional calculus, for studying the behavior of the 
atom-field interaction in photonic crystals. By studying the spontaneous emission of an atom in 
a photonic crystal with one-band isotropic model, we found that the long-time inducing memory 
of the spontaneous emission is a fractional phenomenon. This behavior could be well described by 

■ the fractional calculus. And the results show no steady photon-atom bound state for the atomic 
resonant transition frequency lying in the proximity of allowed band edge which is encountered 
in the previous study [J. Opt. B: Quantum Semiclass. Opt. 5, R43 (2003)]. The correctness of 

■ this result is validated by the "cut-off smoothing" density of photon states (DOS) with fractional 
C ■ calculus. By obtaining a rigorous solution without the multiple-valued problem for the system, we 

m 

show the method of fractional calculus has logically concise property. 
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I. INTRODUCTION 



The performance of photonic devices in various fields is greatly limited by spontaneous 
emission rate. In light-emitting diodes and lasers, for example, spontaneous emission that is 
not extracted from the devices will contribute to loss and noise. Inhibiting undesirable spon- 
taneous light emission and redistributing the energy into useful forms becomes important in 
these fields . It has been demonstrated theoretically and experimentally 0, Q] 

that photonic bandgap (PBG) materials could be effectively used to inhibit the spontaneous 
emission. Near a photonic band edge, the photon density of state (DOS), which determines 
the rate of spontaneous emission, is significantly different from that of free space. Singu- 
larity 81 of the DOS near PBG leads to the strong atom-field interaction and formation of 

n n 

photon-atom bound states [9|, |10[, where the spontaneous emission rate is inhibited. The 
Markov approximation ll| of spontaneous emission in free space is no longer valid in this 
near PBG region, where the atomic decay becomes non-exponential and the emission spec- 
trum becomes non-Lorentzian. By using two-dimensional (2D) photonic crystals, Fujita et 



al. 



12J successfully inhibited and redistributed the spontaneous light emission by a factor 



of 5 clS cL result of the 2D photonic bandgap effect. 

The experimental data clearly indicate that the spontaneous emission rate has non- 
exponential decaying behavior when the emission peak is located near the band edges (lattice 

This non-Markov behavior of the PBG 



131 ] using the Laplace transform method to 



constant a=480 and 390 nm in Fig. 3) of Ref. 
reservoir had been studied by John et al. 
solve the time evolution integral equation of the excited probability amplitude of an atom 
in a high-Q microcavity with singular DOS. 

They showed that the time evolution of the excited-state population exhibits decay and 
oscillatory behavior before reaching a nonzero steady-state value due to photon localization 
@, Q]- This bound dressed state leads to the fractionalized steady-state atomic population 



12 



in the excited state. This behavior is observed as the prolonged lifetime effect in Ref. 
However, John et al. Jij, Q predicted that the unphysical bound state is present even 
when the resonant atomic frequency lies outside the band gap. This is inconsistent with the 
experimental result that the prolonged-lifetime effect will disappear when emission peak lies 
outside the PBG region. This inconsistency may be caused by the multiple- valued problem 
countered in the studies 13] and by the singular DOS which is appropriate for high-Q mi- 
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crocavity but not suitable for PBG reservoir with more smooth DOS near band edge. 

The time evolution of the probability amplitude of excited level of an atom is related to 
the delay Green function or memory kernel G(t — t') ji], 11], which is a measure of the reser- 
voir memory on the excited atom. The resultant Green function depends very strongly on 
the photon density of states of the relevant photon reservoir. The density of states near the 
band edge in the isotropic one-band model has the form of p{oS) oc (uj — co> c ) _1 / 2 , where the 
square- root singularity is a characteristic of a one-dimensional phase space. The correspond- 
ing memory kernel [§| has the same form of square-root singularity: G(t — t') oc (t — t')~ 1 / 2 . 
Such a square-root singularity makes the resultant memory kernel possessing the long-time 
memory effect. That is, there is no time scale to separate the microscopic levels from the 
macroscopic levels. The long-time memory phenomena have attracted a great attention in 
statistical physics. This stochastic property of temporal degree of freedom can be well de- 



scribed by fractional calculus 



lq . Il7j ] . The fractional calculus becomes popular in which 



a diffusion process generated by a fluctuation with no time scale at the macroscopic level. 
This anomalous diffusion process can be further described by fractional Langevin equation 

0,Q. 

In this paper, we applied the fractional calculus to study the dynamics of the spontaneous 
emission of an atom in a photonic crystal. We derived a fractional Langevin equation for this 
system and solved it to obtain the excited-state pro bability density. The solution produced 
by he fractional (inverse) Laplace transform 20|, |2l|, |22j is expressed in terms of the square 
complex variables. There is no multiple- valued problem encountered in the previous studies 



13 



14j. This rigorous mathematical method shows that no steady photon-atom bound state 



exists for the atomic resonant transition frequency lying in the allowed band. We verified the 



correctness of this fact by using the "cut-off smoothing" DOS [23| with fractional calculus 
for the atomic transition frequency lying in the proximity of the allowed band edge. The 
excited-state probability of this result still show decaying characteristics. Fractional calculus 
gives the correct description of the behavior for the system near the band edge either with 
or without the "cut-off smoothing" DOS. It not only resolves the multiple- valued problem 
but also avoids choice of smoothing parameter. Therefore, we suggest that the behavior of 
the atom-field interaction in photonic crystals should be expressed in terms of the fractional 
calculus. 
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II. THE DYNAMICS OF THE SPONTANEOUS EMISSION 



The system we investigate is a two-level atom coupled to the radiation in a photonic crys- 
tal with one-band isotropic model. In a rotating-wave approximation, the total Hamiltonian 
for the coupled atom-field system can be written as 

H = hu 21 a 22 + huj k a t a k + ih J2 9k( a t a ^ ~ ^2iOfc) (1) 
k 

where = 1,2) are the atomic operators for a two-level atom with excited 

state |2), ground state |1), and resonant transition frequency uj 2 i; and at are the annihi- 
lation and creation operators of the radiation field; uj; is the radiation frequency of mode k 
in the reservoir, and the atom-field coupling constant = CJ21 ^ 21 [ 2 2 % • Ud is assumed 

— * 

to be atomic position-independent with the fixed atomic dipole moment d 2 i = d 2 iUd, the 

— * 

sample volume V, polarization unit vector of reservoir mode k, and the Coulomb constant 

co- 
in the single photon sector, the wave function of the system has the form 

W)) = A{t)e-^\2,{Q}) + £3e(*)e-^|l,{l jE }> (2) 

k 

with initial condition A(0) = 1 and B^(0) = 0. Here A(t) labels the probability amplitude 
for the atom in its excited state \2) with an electromagnetic vacuum state and B^(t) for the 
atom in its ground state |1) with a single photon in mode k with frequency u^. 

We got the equations of motion for the amplitudes by projecting the time-dependent 
Schrodinger equation on the one-photon sector of the Hilbert space as 

J t Mt) = -T,9Mt)e'^ t (3) 

jBkit) = 9 % A{ty^ (4) 

with detuning frequency = uj^. — uj 2 i. By substituting the time integration of Eq. (4) into 
Eq. (3), we have the time evolving equation of the excited-state probability amplitude 

jA{t) = - j^G{t-r)A{r)dr (5) 
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with the memory kernel G(t 



T 



-m s (t-r) = p\ J p^e-^-^dco. We could 



observe from this equation that the memory kernel is a measure of the reservoir's memory 
in its previous state. The system evolves proceeding according to this equation. 

Our studying aim, the dynamics of the spontaneous emission, could be obtained by 
solving the time evolving equation (5). Here we applied fractional calculus and fractional 
Laplace transform to solve this equation, which were shown to have the form of fractional 
Langevin equation. We adopted the isotropic one-band model to determine the memory 

o3/2 

kernel in equation (5), which is given as G(t — r) = ( t _^i/ 2 e~^ 7r / 4 ~ Ac ^~ T ^ with t > r in the 
long time limit [8| . Here A c = uj2\ — uj c is the detuning frequency of the atomic resonance 



frequency u 2 i from the band edge u c and P x 



3/2 



7/2 . 2 



21 



d^) / (127r 3 ^ 2 he c 3 ) is the coupling 



constant. By substituting this memory kernel to Eq. (5), making a transformation A(t) 



e iAct C(t), and app 



operator 



21 



ying the definition of the Riemann-Liouville fractional differentiation 
221 ] . we got the fractional form of the time evolving equation as 

1 



J3/2 Jl/2 

' -C(t) + zA c — - 2 C(t) + ^ /2 e-^C(t) 



-f 



-3/2 



(6) 



dt 3 / 2 K> c dty^ K " J ' v "^ i " ~ w 2v^F 
which is defined to be a fractional Langevin equation of this atom-field interaction system. 

We proceeded to solve the probability amplitude C(t) by performing (inverse) Laplace 
transform on this fractional Langevin equation. It gave 



C(t) = £ a n 



71=1 



(?) 



where EAa, a) 



o r(i+n+i) * s ^ ne fractional exponential function of order a, X n are 
the roots of X 3 + iA c X - (i(3f' 2 = 0, and a n are X n -related coefficients. X n and a n are 
expressed as 

X 1 =p 1 /2( v+ + v _) e i«/^ 

X 2 = (3 1 /2( 7]+e -™/z 
X 3 = ^ 2 (rj+e^ 6 - 



r]_e iT/6 )e- i7T/ \ 



r].e-' l7r/6 )e l37T/ \ 



(8) 
(9) 
(10) 



with 



and 



X n 



v± 



1 1 

2 ± 2\ 



v 27 (3 3 ) 



■1/3- 



(X n — Xj)(X n — X m ) 



(n^j^ m;n,j,m = 1,2,3). 



(11) 



(12) 



We have the dynamics of spontaneous emission for the system by plotting the excited-state 
probability density P(t) = \A(t)\ 2 = |C(t)| 2 , which has no multiple- valued problem. This 
probability amplitude could be further expressed by error function as 

A(t) = £ a n {X n + Y n (Erf(^xft)))e x ^ (13) 

n=l 

with Y n = Jx 2 (n = 1,2,3), which is the result of previous studies Q, Q. However, this 
expression has multiple- valued problem because the square roots of complex numbers X 2 will 
introduce multiple- valued complex numbers to Y n . When the numerical results are shown, 
every complex sheet of X 2 has to be checked one by one to avoid errors. Obviously, our 
results of fractional calculus are mathematically rigorous and concise. Besides, our result 
in Fig. 1 shows that there is no steady photon-atom bound state for the atomic transition 
frequency lying in the proximity of allowed band edge which is very different from the results 
of previous studies [3, Q . 

In order to fortify the accuracy of the result of no unphysical photon-atom bound state in 
the allowed band, we proceeded through applying the "cut-off smoothing" density of state 
(DOS) to investigating the behavior of the system near the allowed band edge j^]. As 
mentioned before, DOS in isotropic single band model p(cu) oc (cu — u; c ) _1 / 2 0(u; — co c ) with 
Heaviside step function 6 has a weak singularity for uj — > u c (near band edge) . This singular 
behavior is treate d by "cut-off smoothing DOS" p s (uj) oc lim e ^ fe^~f+^ d(uj — cu c ) in realistic 
photonic crystals 23|. Here e is smoothing parameter and the superscript s denotes the "cut- 
off smoothing DOS" case. The excited-state probability amplitude A s (t) = e lAct C s (t) could 
be exactly solved by performing Laplace transform on the memory kernel and time evolving 
equation (5). These procedures gave the Laplace transform of probability amplitude C s (t) 
as 

° ^ ~ s + iA c + G(s) ^ 
with G(s) = This expression could be further rewritten in terms of the roots, 

X s n {n = 1, 2, 3), of z 3 + Vi~ez 2 + iA c z + (zv^A c - (3 3 / 2 i 3 / 4 ) = as 

where 

I^^^+Oe"' 4 -^ (16) 
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X s 2 = P^ 2 (r] s + e- i77/6 - ri s _e i7T/6 )e~ i7T/4 - ^ (17) 



XI = ^l\rf^l % - nLe-^y**/* - — (18) 



with 



and 



1± ^2 2 J ,X 3 /3 3 / 2 27 /j 3 / 2 ' 4 X 27 V /3 J 

o - ( ^4^^) ( "^ /m; - j -^ 1 ' 2 - 3) (20) 

Here again we used fractional calculus (fractional inverse Laplace transform) to obtain the 
excited-state probability amplitude 

3 



(21) 



A'(t) = e"°< < W-i (K?) + Ke {X - )H 

n=l L Z J 

As we study how the system behaves near the allowed band edge, we choose the detuning 
frequency A c = 0.3/3 with smoothing parameter e = 1CT 3 ,10~ 5 ,0, respectively in Fig. 2. 
It could be observed that the excited-state probability density P s (t) = \A s (t)\ 2 has small 
oscillatory behavior in the short time regime but approaches zero in the long time limit. 
It means that there is really no steady photon-atom bound state for the atomic transition 
frequency lying in the proximity of allowed band edge. Actually we have plotted all the 
behavior of the system inside the allowed band in Fig. 3. For the atomic resonant transition 
frequency located deep inside the allowed band (A c / f3\ = 10) or very close to the band edge 
(Ac//?! = 0.01), we found that these probabilities all show decaying characteristics in the 
long time limit. That is, the photon located within the allowed band will not strongly interact 
with atom so the photon-atom bound state will not be formed. This result could be verified 
analytically from the mathematical expression of the excited-state probability amplitude 
A s (t) = e iAct J2l=i < E t(-\, PQ) 2 ) + X^e (x " )2 ' . For the positive detuning (inside the 
allowed band A c > 0), both terms in the square bracket will asymptotically cancel out each 
other as time approaches infinity {t — > oo). We get correct depiction of the dynamics of 
spontaneous emission in a PBG reservoir through fractional calculus, which is proposed as a 
better mathematical method to study the behavior of the atom-field interaction in photonic 
crystals. 
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III. CONCLUSION 



The dynamics of the spontaneous emission of an atom in a photonic crystal with one-band 
isotropic band structure can be treated by the fractional calculus using either singular or 
"cut-off smoothing" density of photon states. For the first time to our knowledge we show 
that it is a fractal phenomenon that induces the long-time memory of the spontaneous emis- 
sion in the photonic crystal. Solving the time evolving equation of the probability amplitude 
for the system governed by the fractional memory kernel described by singular density of 
states, we obtained rigorous solutions without multiple-valued problem encountered. Be- 
sides, we found that there is no unphysical state of fractionalized atomic population in the 
excited state when the resonant atomic frequency lies in the allowed band, even extremely 
close to the band edge. This result was validated by the "cut-off smoothing DOS" with frac- 
tional calculus. We suggest that the correct description of the dynamics of the spontaneous 
emission in a photonic crystal should be expressed in terms of the fractional calculus. This 
mathematical method, concerned about the isotropic model here, can be easily extended to 
the anisotropic ones and to study interesting effects such as the enhancement of the index of 
refraction with greatly reduced absorption, electromagnetically induced transparency, and 
optical amplification without population inversion. 
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FIG. 1: Excited-state probability density, P(t) = \A(t)\ 2 , as a function of j3t for various values of 
the atomic detuning frequency (A c = U21 — w c ), A c //3 = — 1 (inside the band gap), A c //3 = (at 
the band edge) and A c //3 = +1 (within the allowed band) 
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pt 

FIG. 2: Excited-state probability density, P s (t) = |^4 s (i)| 2 , for the atomic detuning frequency 
A c //3 = 0.3 with three values of smoothing parameter e = (solid line), e = 1CT 5 (dashed 
line),e = 1CT 3 (dot dashed line). The difference of these lines marked by circle is enlarged and 
shown in the inset 
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FIG. 3: Excited-state probability density, P(t) = \A(t)\ , for various values of atomic detuning 
frequency inside the allowed band (A c / /3 > 0) 
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